T.A. Aslanyan 

Institute for Physical Research, Armenian National Academy of Sciences, Ashtarak-2, 378410 Armenia 



in 
o 
o 

(N 
(N 



c3 



13 

o 
o 



> 
oo 

(N 
O 
in 
o 

-I— > 
c3 



i 

O 

o 



X 



_Depi. o/ Appl. Phys. and Chem. 



T. Shigenari and K. Abe 
University of Electro- communications, 
(Received) 



Chofu-shi, Tokyo 182, Japan 



We show that the method of calculations of the Debye- Waller factors for IC structures by J.D. 
Axe, Phys. Rev. (1971), B 21, 4181, is not applicable for such calculations. The set of variables, 
which should be used for such calculations, is the phase ip and the amplitude of the IC modulation, 
but not the real and imaginary parts of the complex IC modulation function. 



The Debye- Waller (D-W) factors for incommensurate 
structures (IC) were first calculated by Overhauseri 
From these calculations followed that the D-W factors 
are tending to zero with approaching the second order 
phase transition, and in some temperature range close to 
the transition the IC satellite reflections should not be 
observed. However, in some tens of the discovered IC 
structures the IC satellites were distinctly observed, and 
their intensity do not tend to zero, as it was predicted by 
Overhauser. The latter became the main motivation for 
the calculations by Axe£, who derived that the D-W fac- 
tors for IC structures are always close to unity, including 
very close ranges near phase transition. However, in the 
paper by Axe^ was not pointed out the reason of such 
a discrepancy between the two results, and we consider 
that this problem is steel in order. We note that the re- 
sult by Axe was reproduced in the book by Bruce and 
Cowleyi, and the result by Overhauser in the book by 
Krivoglaa^. 

Below we compare the calculations of the D-W factors 
by Axe£ with those by Overhauseri and show that the 
variables used by Axe are not applicable for such calcu- 
lations, and therefore the calculations by Axe can not be 
treated as an alternative to those by Overhauser. The 
starting point of the discussion should be the expression 
for the thermodynamic potential: 
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where the condition for the IC transition assumes that 
the coefficients g < and g\ > 0, and the value of the IC 
wavevector k is given by k 2 — —g/2g\. The IC structure 
which appears below the transition may be described us- 
ing two geometrically equivalent forms of the modulation 
function: 



r](R) = (rjo + Srj) cos(kR + (p) 



(2) 



where 770 is the equilibrium IC amplitude, and 5rj is its 
fluctuation. In the equivalent form: 



r)(R) = 771 cos kR 4- 772 sin kR 



(3) 



where 771 and 772 are the real and imaginary components 
of the complex IC wave's amplitude Qk, be., 



and 771 = (Qk+Q-k)/2, 772 = i(Qk.—Q-k)/2, and besides 
Vi +V2 = (Vo + $V) 2 , mhi = tan ip. 

It should be emphasized that for the geometrical de- 
scription of the fluctuation-induced distortions of the IC 
modulation wave the both sets of variables, namely <5?7 
and <p eq.(2) (the IC amplitude and the IC phase), and 
771, 772 eq.(3), are equivalent. However, calculations of 
the statistical averages using these sets of variables give 
significantly different results. The set of the variables Si] 
and ip was used in the calculations of the D-W factors for 
the IC structures by Overhauser ^1 Krivoglaa^ and also by 
the present author A The variables 771 and 772 were intro- 
duced for description of the IC modulation fluctuations 
(phasons) by Bruce and Cowley^ and they were used in 
the calculations of the D-W factors by Axem. In what 
follows we explain why the variables 771 and 772 can not 
be used for such calculations. 

Now we introduce the variables 771 and 772 in the ther- 
modynamic potential for further calculation of the statis- 
tical averages. We should substitute the IC modulation 
function eq.(3) in the thermodynamic potential (1). In 
order to keep our expressions simple we use for the IC 
modulation (3) the form: 



rj(R) = (771 + 5t)i ) cos kR + <5t7 2 sin kR 



(5) 



where 771 is the equilibrium value of 771 and Srji is its fluc- 
tuation. In other words, we assume that the equilibrium 
value of 772 is equal to zero (i.e., 772 = and J772 is its fluc- 
tuation). Therefore, 771 is the equilibrium IC amplitude 
1*71 1 — V01 and the initial phase (po = arctan^/^i) = 0. 
Substituting (5) into eq.(l) one obtains: 
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with<5 = -(a+gk 2 + gi k 4 )/2 > 0, D = {g+Qg 1 k 2 )/2 > 0, 
77J = —a/b. 

The harmonic part of eq.(6) (given by the first line 
of eq.(6)), in fact, was used by Axe- for calculations of 
the D-W factor. As easy to see, this harmonic part is a 
diagonalized quadratic form with respect to Sr/i and ^772 . 
Transferring in eq.(6) to the Fourier-components (using 
only its harmonic part) one can calculate in the harmonic 
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approximation the fluctuations: 

T 1 
V a + Dk z 
T 1 

< mum-k >= (7) 

Since, as it follows from this equations, the mode 772 is 
a gapless mode, it was interpreted by Axe^ as a phason 
mode, and the mode 771 , correspondingly, as an amplitude 
mode. 

Now we demonstrate that the harmonic part of eq.(6) 
can not be used for such calculations in the Landau the- 
ory. Let us calculate the fluctuation correction to the 
harmonic part of eq.(6) from the inharmonic terms of 
the same potential 9 . For example, a correction from the 
inharmonic term bfjiSijiSr]^ can be calculated after pre- 
senting this term in the Fourier components as: 

?7lq?72k?72-k-q 

In the second order of expansion over this inharmonicity 
one will obtain, in particular, the correction 

V 3 f 
A ~ jrb 2 vllmd 2 J dk < \ mk \ 2 >< |77 2k+q | 2 >~ 

The angle brackets <> in eq.(8) denote statistical 
averaging with distribution function exp[— Hq/T], where 
Hq is the harmonic part of eq.(6) in the Fourier-form. 
This correction renormalizes the coefficient at the term 
r\\ q r]\- q in the thermodynamic potential, i.e., for q —> 
it renormalizes the coefficient a in the expansion (6). 
Eq.(8) can be presented also by the graph 9 - 

-q < ^^^ > q 

where the outgoing lines correspond to the mode 771 and 
the internal coupled lines to the mode 772. 

As easy to see, for the wavevector q — > or q = the 
correction A diverges as 

ko 

f dk 1 



regardless of the temperature or closeness to the transi- 
tion. The latter means that the harmonic part of eq.(6) 
can not be used for any calculations, since the inharmonic 
correction to it is not small (it is diverging). One may 
expect that calculation of all the series of the diverging 
graphs would cancel this divergency. However, there is no 
reason to expect that the residual of such a cancellation 
could be neglected. 

We give an example to demonstrate that the calcula- 
tions using variables 771, 772 give physically inconsistent 



results. Since the IC modulation function t](R) itself 
is directly a component of the crystal's electron den- 
sity function it is meaningful to discuss diffraction (of 
X-rays) directly by it, and calculate the corresponding 
D-W factor. Let us calculate the statistical average of 
the modulation function (5) using distribution function 
exp[— Hq/T], where Hq is the harmonic part of eq.(6). 
Since < 6rji >= and < S1J2 >= 0, the result of such av- 
eraging for eq.(5) will be < rj(R) >— 771 coskR. In other 
words, statistical averaging (in the approach adopted by 
Axe) does not reduce (or affect) the amplitude 771 of the 
IC modulation. Since the first order IC satellite reflec- 
tion intensity is proportional to fj 2 , the latter means that 
the D-W factor is always equal to unity for all temper- 
atures. In other words, even at temperatures 10 5 K the 
D-W factor appears to be equal to unity, when it is ex- 
pected to be zero, since at very high temperatures all the 
Bragg-reflection intensity is expected to transfer into the 
diffuse scattering. 

The calculations of the D-W factors by Overhauseri, 
in fact, are based on the averaging of the modulation 
function i](R) given by eq.(2). Taking into account that 
the average^ 

< exp[?:fci? + ip] >= cxp[?:fci?] exp[- < tp 2 > /2] 

one can see that the average of the modulation function 
eq.(2) takes the form 

< r](R) >= 77 coskRcxp[- < ip 2 > /2], 

where the exponent is the D-W factor (it is less than 
unity) which reduces the IC amplitude 770 and subse- 
quently the IC satellite reflections. The local fluctuation 
of the phase < ip 2 > was calculated by Overhauser using 
the thermodynamic potential as a function of the phase 
<p and amplitude 5rj (see also the paper by Golovko and 
Levanyuk4, where a more distinct procedure is used for 
deriving of that potential) . The fluctuation < ip 2 > ac- 
cording to Overhauser diverges as T / rfa with approaching 
the phase transition with 770 — > 0. As a result, the D-W 
factor is tending to zero and attenuating the intensity of 
satellite reflections. The D-W factors for IC structures 
are calculated also in the paper by the present authors 8 
in frame of the Landau theory using the variables ip and 
Srj. Our result significantly differs from that by Over- 
hauser, and it gives not so strong divergence of the local 
fluctuation < tp 2 >~ T 1/2 /r] . 

In conclusion we note that though the use of the vari- 
ables ip and 5rj is consistent with the Landau theory (in 
contrast to the variables Sr/i and £772), the limits of the 
Landau theory applicability to the IC transitions are not 
defined. One may assume that in some region close to 
the transition the Landau theory is not applicable. How- 
ever, there is no reason to expect that the divergence of 
the local fluctuation < ip 2 > will stop on approaching the 
transition. 
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